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Abstract. We give an asymptotic evaluation of the complexity of spherical p-spin spin- 
glass models via random matrix theory. This study enables us to obtain detailed infor- 
mation about the bottom of the energy landscape, including the absolute minimum (the 
ground state), the other local minima, and describe an interesting layered structure of 
the low critical values for the Hamiltonians of these models. We also show that our ap- 
proach allows us to compute the related TAP-complexity and extend the results known 

I in the physics literature. As an independent tool, we prove a LDP for the fc-th largest 

I eigenvalue of the GOE, extending the results of [BDGOl) . 

o 

(N 
> 

O 1. Introduction 

2: 

How many critical values does a typical random Morse function have on a high dimen- 
sional manifold? How many of given index, or below a given level? What is the topology 
,_, of level sets? Our work addresses the first two questions above for a class of natural ran- 

dom Gaussian functions on the iV-dimensional sphere, known as p-spin spherical spin glass 
Oh models in the physics literature and as isotropic models in the Gaussian process literature. 

c~| The third question is covered in the forthcoming paper [AA11| . 

We study here the complexity of these random functions, i.e the exponential behavior 
of the mean number of critical points, and more importantly of the mean number of 
i— —i critical points of given index in a given level set. We introduce a new identity, based 

on the classical Kac-Rice formula, relating random matrix theory and the problem of 
counting these critical values. Using this identity and tools from random matrix theory 
Q\ (mainly large deviation results), we give an asymptotic evaluation of the complexity of 

these spherical spin-glass models. Our study includes the important question of counting 
the mean number of local minima below a given level, and in particular the question of 
finding the ground state energy (the minimal value of the Hamiltonian) . We show that 
this question is directly related to the study of the edge of the spectrum of the Gaussian 
O Orthogonal Ensemble (GOE). 

The question of computing the complexity of mean-field spin glass models has recently 
J> been thoroughly studied in the physics literature (see for example (CLR03] and the refer- 

ences therein), mainly for a different measure of the complexity, i.e. the mean number of 
solutions to the Thouless-Anderson-Palmer equations, or TAP-complexity. Our approach 
to the complexity enables us to recover known results in the physics literature about 
TAP-complexity, to compute the ground state energy (when p is even) , and to describe an 
interesting layered structure of the low energy levels of the Hamiltonians of these models, 
which might prove useful for the study of the metastability of Langevin dynamics for these 
models (in longer time scales than those studied in |BDG0l| ). 

The paper is organized as follows In Sectio n [2} we give our main results. In Section[3j we 



X 



prove two main formulas (Theorem 2.1 and 2.2), relating random matrix theory (specifi 



cally the GOE) and spherical spin glasses. These formulas are consequences of the classical 
Kac-Rice formula (see [AT07J and [AW09J for two excellent recent books giving a very com- 
plete account of this formula and its consequences). The version we use here is proved in 
section 12.1 of the book jAT07] . For another modern account of the Kac-Rice formula and 
similar techniques to those of section [3] see Chapter 6 and section 8.3 of [AW09J). The 
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main ingredient to derive results from the Kac-Rice formula is the fact that, for spheri- 
cal spin- glass models, the Hessian of the Hamiltonian at a critical point, conditioned on 
the value of the Hamiltonian, is a symmetric Gaussian random matrix with independent 
entries (up to symmetry) plus a diagonal matrix. This implies, in particular, that it is 
possible to relate statistics of critical points of index k to statistics of the k + 1-th smallest 
eigenvalue of a matrix sampled from the GOE. 

In Section [4| we compute precise logarithmic estimates of the complexity using the 
known large deviation principle (LDP) for the empirical spectral measure [BG97] and for 
the largest eigenvalue of the GOE [BDGOlJ. In fact we need a simple extension of the 
last LDP, i.e. an LDP for the law of the k-ih largest eigenvalue, which is of independent 
interest and proven in Appendix A. 

In Section [5j we show how these logarithmic results can be used to extract information 
about the lowest lying critical values. We first prove that the lowest lying critical points 



have an interesting layered structure, Theorem 2.15 We then show how our logarithmic 
results imply a lower bound on the ground state energy (the minimal value of the Hamil- 
tonian). At this point it would be useful to have a concentration result for the number of 
local minima, for instance using a control of its second moment. Unfortunately, we cannot 
prove directly such a concentration result. Nevertheless we prove that our lower bound 
is tight (for p even), by proving the corresponding upper bound, using the Parisi formula 
for the free energy at positive temperature, as established by Talagrand [Tal06] . It is re- 
markable that the ground state is indeed correctly predicted by our very naive approach, 
i.e. by the vanishing of the "annealed" complexity of the number of local minima. We 
expect this to be true for all models where Parisi's one-step replica symmetry breaking 
holds at low temperature. In Section [6j we extend our results to the TAP-complexity and 
compare our results to the physics literature [CLR03], |CS95] . 

In Section [7j we show how one can go further and obtain sharper than logarithmic 
asymptotic results for the complexity, using classical tools from orthogonal polynomials 
theory, i.e. Plancherel-Rotach asymptotics for Hermite functions. This section is techni- 
cally involved, so we restrict it to the study of the global complexity, i.e. the mean number 
of critical points below a given level, and we do not push it to include the mean number 
of critical points below a given level with a fixed index. However, we remark that at low 
energy levels, the total number of critical points coincides with the total number of local 
minima. 
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2. Notations and main results 



We first introduce the p-spin spherical spin-glass model. We will fix p an integer larger 
or equal to 2 (the case p = 2 is rather trivial regarding our complexity questions, it will 
be discussed below only in Remark 2.3). 
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A configuration a of the p-spin spherical spin-glass model is a vector of R satisfying 
the spherical constraint 

^S>- 2 = L («) 



N 
i=i 

Thus the state space of the p-spin spherical spin-g lass model is S N ~ 1 {\fN) C R N , the 
Euclidean sphere of radius \fN . 

The Hamiltonian of the model is the random function defined on S N ~ l (^fN) by 

1 N 

H N , p {cr) = N(p _ 1)/2 Y, Jh,-,i P °ii---°i P i <r = K...,(Tjv)e^- 1 (^) 1 (2-2) 
ii,...,ip=l 

where Ji lt ... t i p are independent centered standard Gaussian random variables. 

Equivalently, i?AT,p is the centered Gaussian process on the sphere S N ~ 1 (VN) whose 
covariance is given by 

JV 

E[H NjP (<t)H n , p (*')] = N 1 '? (j2<Ti<rtf = NR(a,a'r, (2.3) 

i=i 

where the normalized inner product R(cr, a') = j^{cr, a') = J2iLi a i a i 1S usually called 
the overlap of the configurations a and a'. 

We now want to introduce the complexity of spherical spin glasses. For any Borel set 
B C R and integer < k < N, we consider the (random) number Crtjv,fc(.B) of critical 
values of the Hamiltonian H^, p in the set NB = {Nx : x G B} with index equal to k, 

Cvt N , k {B)= J2 l{H NiP (a)€NB}l{i(V 2 H N:P (a)) = k}. (2.4) 

(r:VH NtP ((T)=0 



Here V, V 2 are the gradient and the Hessian restricted to S N ~ 1 (i/N), and z(V 2 i?jv,p(f)) 
is the index of V 2 Hn :P at er, that is the number of negative eigenvalues of the Hessian 
V 2 Hn, p - We will also consider the (random) total number Crtjv(-B) of critical values of 
the Hamiltonian Hn, p in the set NB (whatever their index) 

Cvt N (B)= Yl l{H N , p (*)eNB}. (2.5) 

<t:VHjv, p (<t)=0 

Our results will give exact formulas and asymptotic estimates for the mean values 
E(Crtjv,jfc(-B)) and E(Crtjv(-B)), when N — > oo and B, k and p are fixed. In particular, we 
will compute lim logE Crtjv,fc(-B) and lim logECrtjv(-B) as N tends to infinity. 

Before giving the central identity relating the GOE to the complexity of spherical spin- 
glass models, we fix our notations for the GOE. 

The GOE ensemble is a probability measure on the space of real symmetric matrices. 
Namely, it is the probability distribution of the N x N real symmetric random matrix 
M N , whose entries (M^, i < j) are independent centered Gaussian random variables with 



L i-3 

variance 



EM 2 - = (2.6) 
13 2N 

We will denote by Eqoe = Eq Q e the expectation under the GOE ensemble of size N x N. 

Let A^ < A^ < ■ ■ ■ < Aj^_ 1 be the ordered eigenvalues of M N . We will denote by 
Ln = J^j^Q 1 the (random) spectral measure of M N , and by pn(x) the density of 
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the (non-random) probability measure Egoe(-^aO- The function pn(x) is usually called 
the (normalized) one-point correlation function and satisfies 

N-i 



f(x)p N (x)dx 



1 

N' 



-E 



GOE 



i=0 



(2.7) 



We now state our main identity 



Theorem 2.1. The following identity holds for all N, p > 2, k S {0, . . . , N — 1}, and for 
all Borel sets Bci, 



E[CrWB)] = 24 



'-(p-l)TEg 0£ 



-a 



(2.8) 



2(p-l) 

Summing the preceding identities for < A; < iV — 1, we easily find the mean total 
number of critical points given a level of energy and relate it to the one-point function. 

Theorem 2.2. The following identity holds for all N , p > 2, and for all Borel sets BcR, 



E[Crtjv(5)] = 2NJ-{p 



exp 



B 



N(p - 2)x 2 
2~p 



\p N (x)dx. 



(2.9) 



Remark 2.3. We want to discuss here very briefly the trivial case where p = 2. This case is 
easier to understand since there is a simpler connection between the 2-spin spherical model 
and the random matrix theory. In fact, if M is a N x N GOE matrix, the Hamiltonian of 
the 2-spin spherical model can be viewed as the quadratic form defined by the symmetric 
matrix M restricted to the sphere S N ~ 1 (\^N), 

H N)2 (x) = (Mx,x). (2.10) 

Therefore, if p = 2, the 2N critical points of the Hamiltonian are simply the eigenvectors 



of M, while the critical values are the eigenvalues of M. Our formulas (2.8), (2.9) simplify 
greatly for p = 2 and are compatible with this obvious observation. Indeed, for p 
these formulas read 



2, 



and 



E[Crt Ntk (B)]=2P£ OE [\ I i : SB] 
E[Crtjv(B)] = 2Np N {B). 



(2.11) 



(2.12) 



We are now in a position to give our main results about the asymptotic complexity 
of spherical spin-glass models. We will see that the following number is an important 
threshold 

fp — 1 



Let I\ : (—oo 
h(u) 2 



-E„ 



Eoo — E c 
be given by 



>(P) 



P 



(2.13) 



E 2 , 



(z*-El)V*dz 



E 2 



El -log 



E 2 



+ logEoo. 

(2.14) 

Remark 2.4. In [BDGOlJ, it is proved that Ii(u) is the rate function of the LDP for the 
smallest eigenvalue of the GOE with the proper normalization of the variance of the entries 
(more precisely, on its domain h(u) = Ii(—u;E OCJ /2), see (A. 2)). 



We now define the following important functions which will describe the asymptotic 
complexity of the p-spin spherical spin-glass models. 



2 lo g(P 

\ log(p 



ilOg(p- 1), 



p-2 
4(p-l)' 
v- 2 . 

4(p-l)' 



h(u), if u < -Eoo, 

if - Eoo < u < 0, 
if < 



(2.15) 



RANDOM MATRICES AND COMPLEXITY OF SPIN GLASSES 



5 



and, for any integer k > 0, 



1 i og (p - 1) - jg^n 2 — (fc + l)Ji(u), if u < 

if it > 



ilog(p-l)-2= 2 



(2.16) 

2 i( W - ij - , :! " 

We note that p (u), &k,p(u) are non-decreasing, continuous functions on M, with maximal 
values | log(p — 1), | log(p — 1) — respectively (see Figure lk. 

We now give the logarithmic asymptotics of the complexity ofspherical spin glasses. To 
simplify the statement, we fix B = (— oo,u), and we write Crtjv^it) = Crtjv^-B), 

Cit N {u) = Crtjv(B). 



Theorem 2.5. For all p > 2 and k > _^xed, 

1 



lim — logECrtAr ifc (n) = G feiP (ti). 

N—>oo iv 



(2.17) 



0.02 r 



0.01 - 



-0.01 



-0.02 




Figure 1. The functions Qk, P for p = 3 and = (solid), k = 1 (dashed), 
k = 2 (dash-dotted), fc = 10, = 100 (both dotted). All these functions 
agree for u > —E^. 



Remark 2.6. It is straightforward to extend the last theorem to general Borel sets B (see 



Remark 4.1 ). Furthermore, by symmetry, Theorem 2.5 also holds as stated for the random 
variables CrtAr i Ar_;((u, oo)), with I > 1 fixed, if one replaces @k,p(u) by 0/_i iP (— u). 



Remark 2.7. For the local minima, i.e. when k 
rem 



0, the limit formula given by Theo- 
2.5 is precisely the formula given by physicists in [CS95], |CLR 03] . Arguing via a 
TAP approach (to be described below in Section [6]), they derive the following asymptotic 
complexity of local minima, 



where z 



p-i 



9(E) 
E - (E 2 



P 



P 



log 



pz 
~2 



+ 



P 



p2 z 2 



(2.18) 



?(p-l) U/2> 



) 1 / 2 ). In Section 6, we show that, in fact, g(E) 



©o, P (2" 1/2 £). The factor 2" 1 / 2 comes from the fact that in |CS95j the Hamiltonian H 
has a different normalization. 

We also provide an exponential asymptotic for the expected total number of critical 
values below level Nu. 



Theorem 2.8. For all p > 2, 



lim — logE Crtjv(n) = Q p (u). 



(2.19) 
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Remark 2.9. As a simple consequence of Theorems 2.5 and 2.8, one can easily compute 
the logarithmic asymptotics of the mean total number of critical points E(Crtjv(^)) and 
the mean total number of critical points of index k, E(Crtj\rfc(R))- 

Jim ^ logE(Crtiv(M)) = J Iog(p - 1), (2.20) 

lim 1 logE(CrW, fc (R)) = J log(p - 1) - (2.21) 

This agrees with formula (13) of [CS95]. Note the fact that the mean number of critical 
points of index k is independent of k (at least in these logarithmic estimates) 



Remark 2.10. Theorems 2.5 and |2.8| are simple consequences of large deviation properties 
for random matrices, using the results of [BDG01] and of [BG97]. In Appendix |A"| we 
recall the LDP for the empirical spectral measure of the GOE proved in [BG97J, and we 
prove a LDP for the k-th largest eigenvalue of a GOE matrix, extending the results of 
|BDG01| . 

Remark 2.11. The case p = 2 is particular since the total complexity is then always non 
positive. But for any p > 3 the complexity is positive. 

Using our results about complexity, we now want to extract some information about 
the geometry of the bottom of the energy landscape HjVp- For any integer k > 0, we 
introduce = E^ (p) > as the unique solution to (see Figure [T] again) . 

Ok, P {-E k {p)) = 0. (2.22) 

These numbers will be crucial in the description of the ground state and of the low- lying 
critical values of the Hamiltonian Hjsip- It is important to note that, for any fixed p > 3, 
the sequence (-Efc(p))fceN is strictly decreasing, and converges to E OCJ (p) as k — > oo. The 
first result we want to derive is about the ground state energy, which we define as the 
(normalized) minimum of the Hamiltonian Hpf tP 

GS N = ^ T inf Hjs; p (cr). (2.23) 



Theorem 2.12. For every p > 3 



lim inf GS^ > -E {p) 

N—>oc 



Moreover for p > 4 even, 



lim GS" = -E (p) in probability. (2.24) 

Remark 2.13. The lower bound on the Ground State follows from our complexity estimates 
and holds for all p > 3. To obtain a matching upper bound we use the Parisi formula and 
the one step replica symmetry breaking as proved by Talagrand [Tal06j . The Parisi formula 
is proven there for every p but the one step replica symmetry breaking is proven only for 
even p's. It might be noteworthy that if we could go beyond our annealed estimates of 
the complexity we would be in a position to prove directly the one-step replica symmetry 
breaking at zero temperature 



By Theorem 2.12, it is improbable to find a critical value below the level —NEq(p). The 
next interesting phenomenon is the role of the threshold E^ (p) . Namely, it is (even more) 
improbable to find, above the threshold —NE^p), a critical value of the Hamiltonian of 
a fixed index k, when N — > oo. Otherwise said, above the threshold —NE OCl (p), all critical 
values of the Hamiltonian must be of diverging index, with overwhelming probability. 
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Theorem 2.14. Let for an integer k > and e > 0, k(e) be the event "there is a 
critical value of index k of the Hamiltonian H^j )P above the level —N(E 00 (p) — e) ", that is 
BN,k( £ ) = {Grtjv,fc((— -EqoCp) + oo)) > 0}. Then for all k > and e > 0, 

limsup -= logP(Bjv,jfc(e)) < 0. (2.25) 

By the last theorem, all critical values of the Hamiltonian of fixed index (non diverging 
with N) must be found in the band (—NEq(jj), —NE co (p)). We now explain the role of 
the thresholds E k {p). Namely, it is improbable to find critical value of index larger or 
equal to k below the threshold —NE k (p), for any fixed integer k. 

Theorem 2.15. For k > and e > 0, let A^ t f.( £ ) to be the event "there is a critical value 
of the Hamiltonian H^ )P below the level —N{E k {p) + e) and with index larger or equal to 
k", that is A Njk (e) = {^Xfc Crt N<i (-E k (p) - e) > 0}. Then for all k > and e > 0, 

limsup - log¥(A N>k (e)) < 0. (2.26) 



Theorem 2.15 describes an interesting layered structure for the lowest critical val- 
ues of the Hamiltonian Hn jP . It says that the lowest critical values above the ground 
state energy (asymptotically —NEq(p)) are (with an overwhelming probability again) 
only local minima, this being true up to the value —NEx(p), and that in a layer above, 
(—NEi(p), —NE2(p)), one finds only critical values with index (local minima) or saddle 
point with index 1, and above this layer one finds only critical values with index 0,1 or 2, 
etc. This picture was already predicted by physicists for minima |CS95j . [CLR03 and for 
critical points of finite indices |KL96j . In particular, this says that the energy barrier to 
cross when starting from the ground state in order to reach another local minima diverges 
with N, since it is bounded below by the energy difference between an index-one saddle 
point and the ground state, i.e. by N(Eq(p) — E\(pj). 



Remark 2.16. Even though it does not follow immediately from Theorem 2.15 and from 
our results on complexity, it is tempting to conjecture that the minimum possible energy 
of a critical point of index k, normalized by N, should converge to —E k {p) (For k = this 



is the statement of Theorem 2.12 ), while likewise the maximum energy of a critical point of 
index k, once normalized by N, should converge to —Eoc- It is also tempting to conjecture 
that the main contribution to the number of critical points of a finite index k is given by 
those whose energy is asymptotically —NEoo. That is, the number of critical points of 
any finite index with energy strictly below —NEoo should be negligible with respect to 
those with energy near —NEoo (with probability going to one, as N tends to infinity). 
However, near any energy value in E G {—E k {p),—E 00 ) there are still an exponentially 
many critical values of index k. We cannot reach those statements at this point because 
our complexity results concern only the first moment of Crt^ (it) . We would need to 
control the concentration of these random variables. 

In Section [7j we show that the precision of Theorem |2 . 8| can be improved and we derive, 
using asymptotic properties of orthogonal polynomials, the following sharp asymptotics 
of E(Crtjv(w)). 

Theorem 2.17. For p > 3, the following holds as N — > oo: 
(a) For u < —E^ 

ECrt ^) = ( , L IJT^-r 1 i w S ^ 1/2 e Ne ^(l + o(l)), (2.27) 
{2p7ry/ 2 -<j)'(v) + h(v, ^y(v) 



where v = — tt^/^ri) an< ^ the functions h, <p and are given in (7.1), 

dTlOb. 
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(b) For u = -E Q 



-'OO 



ECrtAK-i^) = 2 ^ {0)V ^ N-^e N@ »(- E -\l + o(l)). (2.28) 
3(p - 2) 



(c) ForuG (-^oo,0) 



ECrt^u) = 2 ^ 2p{E ™~ u2) e Ne ^(l + o(l)). (2.29) 
(2 — p)iru 



(d) For u > 



ECrtjv(u) = 2ECrtjv(0)(l + o(l)) = 4 ^ JV 1 / 2 ^ 9 "^^ + o(l)). (2.30) 

yvr(p-2) 

Since G^(n) < 6o(u) for all k > and u < — -E?oo> we obtain as an easy consequence of 
Theorem |2 . 5| the following sharp asymptotics for the mean number of minima. 

Corollary 2.18. For u < -E^, 

(2pvr) 1 /^ ^(m) + ^(t,, ^)'(u) 



where v = ~ u ^j 2{^~Y) an d the functions h, (ft and h(v, are given in (7.1), (7.10). 
3. Proof of the central identity 



In this section we prove Theorems 2.1 and 2.2 For the proofs, we find it more convenient 
to work with processes of variance one on the unit sphere S 1 ^^ 1 C R N rather than to work 
N,p- 



with Hwv Hence, for <x E S N 1 we define 



f N , p (<r) = ~H NjP (VN<t). (3.1) 



We will regularly omit the subscripts N and p to save on notations, / = /jv>- The 
function / is again centered Gaussian process whose covariance satisfies 

N 

E(/(<x)/(<t')) = (j> l0 i) P = i?(<x,<x') P , tryeS N -K (3.2) 

i=i 

To estimate the mean number of critical points in a certain level we will use the Kac- 
Rice formula as it appears in the recent book of Adler and Taylor [AT07j which we now 
formulate as a lemma. We use (x, y) to denote the usual Euclidean scalar product, as 
well as the scalar product on any tangent space T&S N ~ 1 . Let V 2 / be the covariant 
Hessian of / on S N ~ 1 defined, e.g., by V 2 f(X,Y) = XYf - V x Yf. Here V X Y is the 
usual Riemann connection and X, Y 6 TS 1 ^ -1 are tangent vectors. On S"^ -1 we fix an 
arbitrary orthonormal frame field (-E^k^tv, that is a set of iV — 1 vector fields E% on 
S N ~ l such that {Ei(cr)} is an orthonormal basis of T (T S N ~ 1 . We write (\> a for the density 
of the gradient vector {Eif{a))i<i < N and detV 2 /(cr) for the determinant of the matrix 
(V 2 /(^,^i)(o-))i<i,i<Ar- 

Lemma 3.1. Let f be a centered Gaussian field on S N ~ 1 and let A = (U a , ^ a ) a ^j be a 
finite atlas on S N " 1 . Set f a = f o VP" 1 : ^ a {U a ) C M 7 ^ 1 -> R and define /f = df a /dxi, 
fij = d 2 f a /dxidxj . Assume that for all a G I and all x,y G f„(i7 a ) the joint distribution 
°f (ff( x )> ftj( x ))i<i<j<N is non- degenerate, and 



1 ' 1 ' l x 



max Var(/S(a:)) + Var(/£(y)) - 2 Cov(/?(x), /£(*/)) < K a \ In \x - y\\- L -P (3.3) 
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for some (3 > and K a > 0. For a Borel set B CM, let 

Cvt f N k (B) = ^(^/(o-)) = k, f(a) G B}. (3.4) 

<t:V/(<t)=0 

Then, using dcr to denote the usual surface measure on S 1 ^^ 1 , 

ECvt f Nk (B) = [ E[|detV 2 /(<r)|l{/(<r) E B,i(V 2 f(a)) = k} \ V/(<x) = 0]<£ CT (0)d<r. 

(3.5) 

Proof. Assumptions of the lemma, which are taken from Corollaries 11.3.2 and 11.3.5 of 
[AT07J, assure that / is a.s. a Morse function and its gradient and Hessian exist in I? 
sense. The lemma can be then proved using the same procedure as Theorem 12.4.1 of 
|AT07| . Our formula (3.5) is analogous to the display just following formula (12.4.4) of 
|AT07j . modulo the term (— l) fc which is missing in our settings since we are interested in 
the number of critical points of / in B and not in the Euler characteristic of the excursion 
set. □ 



An application of Lemma 3.1 is made possible due to the following lemma which de- 
scribes the joint law of Gaussian vector (/(<x), V/(er), V 2 /(cr))). A similar computation 
is also present in section 8.2 of |AW09j . 

Lemma 3.2. (a) Let (fiicr))^^]^ be the gradient and {fij{cr))i<i,j<N the Hessian matrix 
at a G S N ~ l , that is ft = ^/(tr),/^- = V 2 f{E i ,E j ){cr). Then, for all 1 < i,j,k < N, 
f(cr), fi{<j), fjk(c) are centered Gaussian random variables whose joint distribution is 
determined by 

E[/H 2 ] = 1, = nfiWfjki*)} = o, 

E[/(<r)ji»] = -p6 ij} mWM")] = rfij, [ ' j 

and 

mij(°~)fki(v)] = p{p - l)(o~ikO~ji + SuSjk) +p 2 5ij5 k i. (3.7) 
(b) Under the conditional distribution P[-|/(er) = x], x G M, the random variables 
fij(cr), 1 < i,j < N, are independent Gaussian variables satisfying 

E LAi(o")] = -xpSij, 
E[^(<x) 2 ] ={1 + Sij)p(p-1). 
Alternatively, the random matrix (fij{cr)) has the same distribution as 

M N ~ l ^2{N - l)p{p- 1) - xpl, (3.9) 
where M 1 is (N — 1) x (N — 1) GOE matrix given by (2.6) and I is the identity matrix. 

Proof. Without loss of generality we can suppose that <x is the north pole of the sphere n = 
(0, ...,0,1). We define the function f : S^" 1 ->• K^" 1 by *(xi, . . . , x N ) = (x x , . . . , £jv-i). 
It is a chart in some neighborhood U of n. We set 

f = foV-\ (3.10) 

which is a Gaussian process on ^(U) with covariance 

,N-l , P 

C(x,y) = Cov(/»,/(y)) = £ x iVi + yj (l - Efi 1 ^)^ - Y.ti vf) ■ (3.11) 

i=l ' 

We choose the orthonormal frame field (E{) such that it satisfies Ei{n) = d/dxi with 
respect to the chart ^> . Then the covariant Hessian (fij(n)) agrees with the usual Hessian 
of / at 0, by noting that the Christoffel symbols rj^(n) = 0. Hence, to check (3.6), (3.7), 
we should prove analogous identities for /(0), fi(0) = gfr/(0) and /y(0) = dx d dx . f(0)- 
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The covariances of /, /j, fij, can be computed using a well-known formula (see e.g. |AT07] 
formula (5.5.4)), 



Cov 



d k f(x) d e f(y) 
dx h . . . dx ik dy h . . . dy je 



d k+e C(x,y) 



dx ix ...dx ih dy h ...dy 



(3.12) 



j i 



Straightforward algebra then gives (3.6), ( |3.7| ). Moreover, since the derivatives of a cen- 
tered Gaussian field have centered Gaussian distribution, relations (3.6) and (3.7) deter- 
mine uniquely the joint distribution of f(cr), fi(cr), and fij(cr). This completes the proof 
of the claim (a). 

The well-known rules how Gaussian distributions transform under conditioning (see, 
e.g., [AT07], pages 10-11) then yield the claim (b) of the lemma. □ 



The next tool in order to prove Theorem |2.1| is the following lemma which is an inde- 
pendent fact about the distribution of the eigenvalues of the GOE. It allows us to deal 
with the (usually rather unpleasant) absolute value of the determinant of the Hessian that 
appears in ( |3.5| ). 

Lemma 3.3. Let M w_1 be a (N — 1) x (N — 1) GOE matrix and X be an independent 
Gaussian random variable with mean m and variance t 2 . Then, for any Borel set GcM, 

rN ~ x - Xl\l{i{M N - x 



E 



| det M 1 

r(f)(iv 



XI) = k,X G G} 



1 



' 7T 



t 2 



¥ N 
GOi? 



exp 



(3.13) 



Proof. The left-hand side of (3.13) can be written as 

i / (x—m) 



G 



N-l 
GOE 



det(M - xI)\l{i(M - xl) = k) 



dx. 



(3.14) 



Observe that the event {i(M — xl) = k} is equal to the event {A^ (%)}, where A^ (x) is 
given by 



{A"" 1 : An -1 <•••! A,: ",' < ,r< A; 1 "' < 



.N-l 



<N-1 



.N-l 



. < A 



jv-r 



, N _fi. (3.15) 

Recall the explicit formula for the distribution Qn of the eigenvalues A^ of the GOE 
matrix M N , (see, e.g., |Meh91| . p. 519), 

N-l 

Q N (dX N ) = ^He-f^ 2 dXf J] \^-^\H X o<--<^-i}, (3-16) 

N i=0 0<i<j<N 

where the normalization Zjv can be computed from Selberg's integral (cf. .Me hOI j . p. 529) 



i N 
Z N = l.^V2) N N- N( - N+1 ^ [J F ( 1 + 

i=l 



Using this notation, 



E 



N-l 
GOE 



det(M - xI)\l{i(M - xl) = k) 



N-2 

TT \Xf~ l - x\Q N -x(dX 

W '=0 



N-U 



and (3.14) becomes 
1 



V2vrt 2 



(x — m)^ 



G 



N-2 

n i^ -1 



x\Q N -x(dX N - 1 )dx 



(3.17) 



(3.18) 



(3.19) 



Comparing the product in the integrand with the van der Monde determinant in (3.16) 
suggests considering x as the fc+l-th smallest eigenvalue of a N x N GOE matrix. Indeed, 
if we substitute A 



N-l 



{^} 1/2 Af, for i € {(),..., A: -1}, A 



N-l 



{^} 1/2 Af +1 for 
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i G {k, . . . , N — 2} and x = {^rj} 1//2 A^ and perform the change of variables, then we 
obtain 



Zjv_ iv^rt^ ViV- 1 



(JV+2)(iV+l) 
4 



x E 



AT 

GOE 



exp 



(^_«^M) lKe{V} v M ], 



Plugging (3.17) into (3.20) completes the proof of the lemma. 



(3.20) 



□ 



The three above lemmas yield the proof of Theorem 2.1 as follows. 



Proof of Theorem 2.1, We first verify the conditions of Lemma 3.1 We use the same chart 
and orthogonal frame as in the proof of Lemma 3.2 From (3.6), (3.7), it is not difficult 
to check that the joint distribution of (fi(cr), fij(cr)) is non-degenerate for a = n. By 
the continuity of the covariances, it is then non-degenerate for all <x G U, if U is small 
enough. Similarly, using (3.12) we can verify that (|3.3l) is satisfied on U. Since S N ~ l can 



be covered by a finite number of copies of U obtained by rotations around the center of 
sphere, the conditions of Lemma 3.1 are satisfied. 

We can thus apply formula (3.5 ). First note that, due to the rotational symmetry again, 
the integrand does not depend on a. Hence, recalling that n denotes the north pole and 

(3.21) 



using (2.5),(3.1), 



ECrtjv,fc(£) 



det V7(n) l{i(V z f N (*)) = fc}l{/(n) G VNB} V/(n) = 0„(O) 



where lon is the volume of the sphere S 



N-l 



UJ N 



T(N/2) 



(3.22) 



The density of gradient V/(n) is same as the density of (/i(0))i<i<jv- Hence, using (3.6), 

MO) = (2vrp)-( JV - 1 )/ 2 . (3.23) 



To compute the expectation in (3.21), we condition on f(n) and use the fact that, by 
(3.6), both / and its Hessian are independent of the gradient, 

\detV 2 f N (n)\l{i(V 2 f N (a)) = kJ(n) G VNB}\Vf(n) = 



E 



E 



(3.24) 



E[\detV 2 f N (n)\l{i(V 2 f N (a)) = k,f(n) G VNB}\f(n) 

By Lemma 3.2, the interior expectation satisfies 

E[| detV 2 / 7V (n)|l{i(V 2 / 7V (cr)) = A;,/(n) G VNB}\f(n)] 

det (M N ' 1 - p l ' 2 {2(N - l)(p - l))-^ 2 f(n)l) \ 



(2(iV_i) p ( p _i ) )^ 1 E^ 



x 1 



{i^" 1 - p 1 / 2 (2(iV - l)(p - l))- 1 / 2 /^)/) = k, f(n) G VNB] 



Inserting (3.25) into (3.24), we can apply Lemma 



3.3 



with m = 0, t 2 



(3.25) 

2(JV-l)(p-l)' an< ^ 

2(jv-i)(p-i) Using (3.22) and (3.23), we get after a little straightforward algebra, 



G — 1 ^ 



2^(p-l)T E ^ OE 



UK g 



2(p- 



tt4 



This completes the proof of Theorem 2.1 



(3.26) 
□ 
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Proof of Theorem 2.2. Theorem 2.2 follows from Theorem |2.1| by summing over k £ 



{0, . . . , N—l}. The additional TV" in the prefactor comes from the fact that pn is normalized 
one-point correlation function, cf. (2.7). □ 



4. Logarithmic asymptotics of the complexity 



In this section we apply the LDP for the k-th largest eigenvalue of GOE, Theorem |A,1[ 
to study the logarithmic asymptotics of the complexity of the spherical spin glass, that is 



we show Theorems 2.5 and 
Theorem 



A.l 



2.8 Observe, that comparing (2.6) with (A.l), we must use 
with a = 2 -1 / 2 . Here and later we write A; for A; . 



4.1. Proof of Theorem 2.5 We analyse the exact formula given in Theorem 2.1 



or 



equivalently in (3.26). By Theorem A.l and the obvious symmetry between the largest 



and the smallest eigenvalues, the (k + l)-th smallest eigenvalue of M N satisfies the 
LDP with the good rate function Jk(u) = (k + u; 2~ 1 / 2 ), wh ere I\ is defined in 



(A. 2). Set t = uw 2 (p-i) an d 



lim — log! Crt jv,fc(ii) = - log(p 



■^2f^x 2 . Then, by Theorem |2.1 

1)+ lim ^logE& E [e*«l Afc < t ] 

iv— »oo iv 



(4.1) 



By Varadhan's Lemma (see e.g. | DZ98| , Theorem 4.3.1 and Exercise 4.3.11, observe that 
4> is bounded from above, so that condition (4.3.2) of [DZ98J is obviously satisfied) 



sup ( 

x&{— oo,t) 



J k {x)) < hminf llogE^ OE [e^)l A; 

A<— >oo iv 



< lim sup i logEgo E [e^ A *)l Ajk < t ] < sup {<f>(x) - J k (x)). 

N— >OC I* ^Ct—rr, fl 



(4.2) 



Xd ( — oo,t] 



It can be seen easily that for t < —y/2 the both suprema in (4.2) equal (f>(t) — Jk{t)- On 



the other hand, if t > — v2, these suprema equal (j)(y/2). Hence, using the definitions of 
4>, t, Jfc, Eoo, and the scaling relation (A.5), 



(4.1) 



\ log(p 
5 log(p 



p-2 

P ' n 

(p~2)u 2 
4(p-l) 



(k + l)h{-u; En/2), 



if u > -Eoo, 

if U < -Eon. 



(4.3) 



Using Remark 2.4, this completes the proof of Theorem |2. 5 



Remark 4.1. The proof of Theorem 2.5 clearly extends to a Borel set B. In fact, one just 
need to apply Varadhan's Lemma and find the supremum of <p(x) — Jfe(x) on a appropriate 
domain. 



4.2. Proof of Theorem 2.8, Let t and <ft as in the previous proof. By Theorem 2.2 
have to study 



we 



1 12 
lim — log2N\ -(p 

N^ooN & VP 



1 TE 



> N ^L N (dx). 



(4.4) 



where Ljq is the empirical spectral measure of the GOE matrix (see below ( |2.6| )). The 
constant in front the integral gives the term 2log(p — 1) of @ p (u). We need to evaluate 
the contribution of the integral. For t < —y/2, using l Ao <t > lAi<* for all i, we write 



^[e^'l Ao < ( ] 



< E 



, N ^L N (dx)<E[e N ^h Xo < t }. 



(4.5) 



Taking the logarithm and dividing by N, the both sides of this inequality converge to 
(f>(t) — t; 2 -1 / 2 ), by the same argument as in the previous proof. Using the values of 
t and (j), this proves the theorem for u < —-Eoo- 
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For t G (-V2,0], we cannot use the smallest eigenvalue Ao in the lower bound. Therefore 
we write, for e > 0, 

e N ^L N (dx) <E[e^Wl Ao < t ]. (4.6) 



1 

N 



E[e«l t , ((1 . £il)) >„] 



< E 



By the LDP for Ao, P[Ao >t] — >• 1 as N — > oo. Similarly, by the convergence of L n to the 
semi-circle distribution, we have P[Ljv((t — e,t)) > 0] — > 1 as N — > oo. Therefore, after 
taking the logarithm and dividing by N, we find 



4>(t-e) < lim ^logE 



i N ^L N (dx) < <f>(t). 



(4.7) 



Since 4> is continuous, the claim of the theorem follows for t G (— y/2, 0], or equivalently 
for u G (— Sjxj, 0]. The proof in the case u > is analogous and we left it to the reader. 

5. The geometry of the bottom of the energy landscape 

In this section, we use our complexity estimates to obtain information about the bottom 
of the energy landscape. We first prove Theorems 2.14 and 2.15 We then prove Theorem 

malized energy of the groi 

1 m£{H Nj> (a) : a G S N -\VN)} (5.1) 
converges to —Eq(p) in probability. 

5.1. Proof of Theorem |2.14| , We want to show that there are no critical points of a 
finite index of the Hamiltonian above the level N^—Eqo + e). Let k and e be as in the 



2.12, that is we show that the normalized energy of the ground state 
GS N = N- 



statement of the theorem. Then, by Markov's inequality and Theorem 2.1 



P[Crtjv,fc((-£;oo + e,oo)) > 0] < E[Crt 7Vife ((- J B 00 + e,oo))] 

iV(p-2)(A fc ) 2 ' 



< c(p)(p- 1) N/2 E 



exp 



2p 



-}l{A fe > -V2 + e'} 



(5.2) 



< c(p)(p-l) N/2 F\\ k > -V2 + e' 



where e' = eyf^jE^. By the LDP for the empirical spectral measure Ln (see Theorem 
1.1 of |BG97j ). we know that for some C(e') > 

-CrN 2 



< e" 



(5.3) 



Combining this estimate with (5.2) completes the proof of Theorem 2.14 



5.2. Proof of Theorem 2.15| , We want to prove that there are no critical values of index 
k of the Hamiltonian below —N(Ef : + e). Using Theorem 2.5, we have 

E[Crt N , k (-E k - e)] <exp{NQ k)P (-E k -e) + o(N)}. (5.4) 

The function Q ktP is strictly increasing on (— oo, —Eoq). The constant E k > E^ is defined 
by @ k ,p(—E k ) = 0. Therefore, ® k:P (—E k — e) = c(k,p,e) < 0. An application of Markov's 
inequality as before completes the proof of Theorem |2.15| 

5.3. Proof of Theorem |2.12[ Note that, by Theorem |2.15 there are no minima of the 
Hamiltonian below —N(Eq(p) + e), with high probability. This implies that for all e > 

lim ¥[GS N > -E (p) - e] = 1. (5.5) 



Y^oo 



To find a matching upper bound on GS N we use known results about the free energy 
at positive temperature, more precisely the Parisi formula as proved by Talagrand [Tal06j. 
Recall that the partition function of the p-spin spin glass is given by 



ZnM 



s n ~ 1 {Vn) 



(5.6) 
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where A^r is the normalized surface measure on the sphere S N 1 (\HV). By Theorem 1.1 
of |Tal06| . 

lim ±-ElogZ N , p (/3)=F p (/3), (5.7) 

iV— »oo iv 

where F p (f3) is given by the variational principle as in (1.11) of |Tal06] . Furthermore, it 
is known that in the case of the spherical p-spin model, the F p {j3) can be computed using 
the following simpler variational problem (see (1.16) in |PT07| ). 

Fp{P)= mf , \ {/3 2 (l + (m- l)/) + (l-I) log(l -q) + - log (l-g(l-m)) }. (5.8) 

q,m£[0,l\ Z \ 771/ 771 J 

We now analyze this variational problem. We replace q = 1 — dj3~ l and m = in 
the last equation, and we denote the function inside the infimum by P(/3, c, d), 

P(P,c,d) = \{(3 2 {l + (c/T 1 - 1)(1 - d^y) 

2 1 , (5-9) 

+ (1 - c" 1 /?) log (1 - (1 - dp- 1 )) + c-^log [1 - (1 - dp- 1 )^ - c/T 1 )] }• 

The following lemma shows that this is the right scaling as f3 tends to infinity. 

Lemma 5.1. There exist constants < e < M < oo, such that, for all (3 large enough, 

FJ/3)= inf P(p,c,d)= inf P(p,c,d). (5.10) 

c,de[0,/3] c,de[e,M] 

Proof. It follows directly from Lemma 3 in |PT07| . □ 

As f3 — > oo, the function f3~ l P(f3,c,d) converges uniformly on the compact set c, d G 
[e, M\. Therefore the last lemma implies that 

lim /3" 1 F„(/3) = inf -\c + pd+ -(log(c + d) - logd)} =: inf P(c,d). (5.11) 

J->oo c,d£[e,M] 2 1 C > c,de[e,M] 





Lemma 5.2. We have that 



7= inf P(c,d) = E (p). (5.12) 

c,de[e,M] 



Proof. The constant Eq(p) is the unique solution of the equation &o tP (— x) = (see (2.16), 
(2.22)). Therefore, to prove the lemma it suffices to show 

e 0lP (-7) = 0. (5.13) 

A critical point of P(c, d) satisfies 

p- 1 = d(c + d) (5.14) 

c 2 {c + d) + c= (c + d)(log(c + d) - log(d)). (5.15) 

Writing y = c + d, it follows from (5.14) that d = (py) -1 , c = y — d = (py 2 — l)/py. From 
(5.14) and d > e, we further obtain d = |(— c + \f c 2 + 4p _1 ). Therefore c,d > e implies 

y = c + d>p- 1 /2_ ( 516 ) 



Inserting these computations into (5.15), we obtain that y is a solution of 



/py -iy py -1 2 _ 1/2 

M/H = ylog(py), y>p 7 . (5.17) 

V vv / vv 



'py 1 - 1\2 py^ - 1 

y + 

py J py 

Setting a = py 2 , this is equivalent to 

g p (a) := (a - l) 2 + p(a - 1) -paloga = 0, a > 1. (5.18) 
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The function g satisfies g p (l) = g' p (l) = 0, g p {l) < and g" is increasing on [l,oo). 
Therefore, g p has a unique minimum ao on [1, oo) and is strictly increasing on [ao,oo). 
Moreover, since g p (p — 1) < 0, the equation ( |5.17| ) has a unique solution satisfying 

(5.19) 



Using the definition of 7 and (5.17) 

/-I 1 

7 = " I — + " + 



py 2 > p — 1 and thus y > \f (p — 1) /p. 

log(py ) 



l/py_ 
2 V 



1 



2/ + 



.'//' 



(5.20) 



To compute Oo, p (— 7) observe that, by (5.20), 7 — E t 



(y—^p-^) 2 - Hence, (5.19) implies 



7 



+ \R 



El = 2y. 



Inserting these results into definition (2.16) of Bo, p and using equation (5.17) yields after 
a little algebra that 0o,p(— 7) = 0. This completes the proof of Lemma 



5.2 



(5.21) 

er 
□ 



We can now prove the upper bound on GS . Note that 



1 



E log Z 



1 



fiN 

Taking the limits N 



N 



-E log 



/3N 

00 and then (3 



f _ e-^^A^der) < -EGS N . (5.22) 

Js n (Vn) 



s n (Vn) 

00, using (5.11) and Lemma 5.2, we obtain 



EGS N < -E (p) + e for N large enough. 



By Borell-TIS inequality (see Theorem 2.7 in [AW09| L 

PflGS^ + EGS^I > e] < e 



Ne 2 



Combining (5.23) and (5.24), we get that for all e > 0, for N large enough, 

F[GS N < -E {p) + 2e] > 1 - e. 
This combined with the lower bound (5.5) completes the proof of Theorem 2.12 



(5.23) 
(5.24) 
(5.25) 



6. The Thouless-Anderson-Palmer complexity 

In this section, we study the mean number of solutions of the Thouless-Anderson- 
Palmer (TAP) equation (6.4), called the TAP complexity. Using the previous results 
of this paper, we give, in Theorem 6.1, a formula for the TAP complexity at any finite 
temperature. In physics literature [CS95j, the TAP complexity was predicted and used to 
derive a formula for the complexity of the spherical p-spin model. In Lemma 6.3, we show 
that the formula of [CS95! agrees with our Theorem 2.5 For further physics interpretation 
of the TAP solutions, such as connections to metastable states, the reader is invited to 
check |KPV931 IC595J and the references therein. 

Let B(0, v^V) C M> N be the open ball of radius centered at 0. We define the 
TAP functional as Kurchan, Parisi and Virasoro in [KPV93J. For m = (mi, . . . ,171^) in 
5(0, y/N) andq=± m 2 £ [0, 1), let 



F TA p(m) 



1 



N 

2 l/2 iV (P+l)/2 ^2 

il,...,i p =l 



Jh.-ipini! ■■■m ip + B Pi p(q), 



where, as usual, Ji lt ...i p are independent standard normal random variables and 



B(q) = B P}0 (q) 



-i]o g (l- 9 )-^(l + (p-l)/ 



pq 



(6.1) 



(6.2) 



Notice that the TAP functional -Ftap can be written in spherical coordinates: Defining 



q-^m G S N - 1 (y/N) and h N>t 



F TAP (m) = fcAp(q,tr) = 2- l / 2 qPl 2 h N , p {a) + B pS {q). 



(6.3) 



RANDOM MATRICES AND COMPLEXITY OF SPIN GLASSES 



16 



The TAP equations are equations for critical points of the TAP functional, 

d 

■F TA p(m) = 0. 



dm. L 

Since q = is not a critical point of -FtaPi the equations (16.41) are equivalent with 



_d_ 
da, 



f TAP (q,a) = 



_d_ 
da,, 



0_ 
dq 



/TAP(<Z, <t) = 2- 3 / 2 pqP/ 2 - 1 h NtP (a) + B'(q) = 0. 



(6.4) 

(6.5) 
(6.6) 



Therefore, a solution of the TAP equation (6.4) must be a critical point of the Hamiltonian 
H]\r p such that q satisfies (6.6). A critical point (q, a) of the TAP functional of index k is 



thus either a critical point of H^^ of index k satisfying (6.6) and ^/tap(q, 0") > 0, or a 
critical point of H^^ of index k — 1 satisfying (6.6) and ^ ,/tap(3S f ) < 0. 



Let Mk(u, f3) represent the number of critical points of index k of the TAP functional 
with normalized energy tiN tP smaller than u at temperature /3" 1 , 

Af k (u,P)= l{h N , p (a)e(-oo,u},i(V 2 f TAP (q,*)) = 0}. 

(9,o-):V/tap(<?,o-)=0 

For each u 6 (— oo, — E^], p > 3, and (3 > 0, we further define 



(6.7) 



(3(u) 



P 



P 



£-2 

2 



U 



2 3/2( ;p _ 1 )Vp-2, 

u*(/3) = sup {v e (-co, -Eoo] : f3(v) < /?}. 
Observe that (5{-) is an increasing function with lim M _ s ._ 00 f3(u) = 0. 

Theorem 6.1. For all (3 > and p > 3, 

(a) Ifu < -E^, k = 

lim i logEAA (n,/3) = 6 0j ,(«*(/3) A «). 

N— >oo iV 



(6.8) 
(6.9) 



(b) IfuK-E^, k > 0, 



1 



lim -logEA4(«,/3) = 8 fc _ ljP («*(/3) At*). 

iv— s>oo jv 



.10) 



(6.11) 



(c) J//3 < ^(--E(fe_i)vo(p)) or if u < --E^-LjvoO 5 ); then M k (u, f3) tends to zero in 
probability as N —> oo. 

Proof. We start by proving parts (a) and (b). Computing -B'(g) and multiplying the both 
sides of ( |6.6| ) by 2 3 / 2 /?(l — g)/p, (6.6) is equivalent to 

2 i/2 



2~ 1 / 2 (p - 1)/3 2 ((1 - g)^ 2 - 1 ) 2 + h N>p (<r)P(l - q)q^ 2 - 1 + 



0. 



P 



Setting 
we obtain 



z = z(q) = (1 



tP/2-1 



1 



(6.12) 

(6.13) 
(6.14) 



2 i/2(p_i)V -w/- V'W 
Thus, a solution (<r, g) to the TAP equation ( |6.4[ ) is a critical point <x of the Hamiltonian 
that satisfies (6.13) and (6.14). 



The next lemma counts the number of solutions of equation (6.13). 
Lemma 6.2. The function f(q) = (1 — q)q p / 2 ~ 1 — A, has exactly two, one or no zeros 



in 



q e [0, 1) ifO < A < \{?f)^, A = y*ff? and A > |(^)^, respectively. 
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Proof. If A = 0, then clearly and 1 are the only zeros. Since f'(q) = q p / 2 ~ 2 (p—pq — 2))/2, 
there is a unique critical point of / in (0, 1), a maxima at q = p ^-. Now, varying A gives 
us the result of the lemma. □ 



Applying the lemma to equations (6.13), (6.14), we see that in order to (cr,q) be a 
critical point, hN tP (cr) must satisfy (ignoring the equalities which have zero probability) 



< 



2 1 /2 / 3( p _i) 



h N , p (<r)±Jh%J*)-El) < 



2(p-2 



p \ p 



p-2 



which is equivalent to (cf. (|6.8[), (6.9)) 



h N>p (tr) < u*(P). 



.15) 



.16) 



Furthermore, if (6.16) is satisfied, the equations (6.13), (6.14) have two solutions q\, q2 
such that < q\ < (p — 2)/p < qi < 1. 



At a critical point, using (6.12) to compute fo/v>, and (6.13) to simplify the terms 
containing q p , 

d 2 

q^2 f TAP (q,cr) 

= P(P A^-2u , 1 f3q p ~ Z p{p - 1) , 



2 V2 

pjp ~ !) 
8 g /3(l - q) 2 



1 )q 2 h NjP ((T) + 



2/3(1 ~qf 



-((p-l) 9 -p + 2) (6.17) 



P 



P 







"2 2 



PiP ~ 1) 



Since, z{q\) = zfa), the above expression is positive for one of qi, q2 and negative for the 
remaining one. Hence, every critical point er of index k of the Hamiltonian with h^^cr) 
as in (6.16) contributes one critical point of /tap with index k and one of index k + 1. 
Therefore, for every k > 0, /3 > and u < —E^, P-a.s., 

A4(n,/3) 



G [1,2], 



.18) 



CrtAr ifc _i(n A u*(/3)) + Crt fc (u A u*(/3)) 

where we defined Crt7v i _i(— oo, c) = 0. Claims (a) and (b) then follows directly from 
Theorem 2.5 using the observation ®k,pW) > 8fc+i )P (tt) for all k > and u < —E^. 
Claim (c) of the theorem then follows from (a), (b) using the Markov inequality. □ 

6.1. Complexity of the minima: derivation of formula (11) of [CS95j. In this short 
section, we verify that the results obtained in Theorem 2.5 for minima, resp. maxima, agree 
with the formula proposed by A. Crisanti and H.-J. Sommers in [CS95J. 



Lemma 6.3. Let z 

Proof. First, note that for u 



M~ u - \A 2 - ^F) and u * - v^- 



then 



1 (2 
2 



V 



i ,P Z \ P — 



p 2 z 2 



(6.19) 



2(P-1) 



, we have that, by (2.16), 



1 



9o,p(v / 2u) = -log(p- 



P 



(6.20) 



which agrees with the left-hand side of (6.19), taking z 



p(p-i) • 



Second, the derivative with respect to u of the left-hand side of ( |6.19[ ) is given by 

@' , p (V2u) = {p- 1)" 1 {(2 - p)u + y/p(pu*-2p + 2)}. (6.21) 
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On the other hand, the derivative of the right-hand side of (6.19) equals 



\ ((4(p - 1) 2 (-1 - u/^-(2(p-l))/p))/(p 2 (-u - vV-(2(p-l))/p) 3 ) 
- (2(-l - u/y/u*-(2(p-l))/p))/(-u - jv? - (2(p - l))/p) 
+ ((-1 - u /^ u 2-(2(p-l))/p)(-u - V« 2 -(2(p-l))/p))/(p - 1)) . 



.22) 



After a lengthy, but straightforward, simplification, the last expression coincides with 



0' o (y/2u). This completes the proof. 



□ 



7. Sharper estimates of the mean number of critical points 
In this section we prove the sharp asymptotic results on the complexity of spherical p- 



spin spin glass, that is Theorem 2.17 To this end we analyse the exact formula ( 2.9 ) proved 



in Theorem 2.2 For u > —E^, a simple application of Laplace's method is sufficient. 
For u < —Esq, more care is needed. We write the one-point correlation function pn{x) 



appearing in (2.9) as a function of the Hermite polynomials and to use the asymptotics of 



these polynomials given by the Plancherel-Rotach formula. 

7.1. Sharp asymptotics in the bulk. We start with the case u > —E^ that is parts 

(p - 2)x 2 



(c), (d) of Theorem 2.17 Setting 



C NjP = 2NJ*(p-l) N / 2 



P 



2(p-l)' 



<j>{x) 



2p 



(7.1) 



we rewrite the expectation of the global complexity (2.9) as 



ECrt N {u) = C N , P e N ^ x) p N (x)dx = C N , p e N<f>{x) p N {x) dx. (7.2) 



For the last equality we used the fact that <j> and pm are even functions. The case u > — Eqq 
corresponds to v < y/2. Since pn converges uniformly to the density of the semi-circle 
law p(x) = 7r~ 1 \/2 — x 2 l| x | <v ^, by Laplace's method, the principal contribution to the 



integral (7.2) comes from the boundary point x = v. Since F'(v) = iff v = 0, we have 
' "\ ,< .»('>), ■— ( 1 + o( I)). 

(7.3) 



ECrt7v(0) = ^C NtP e mo) p(0) 



2 ~"*~ — y NFg(Qy 

ECvt N (u) = C N , p e N ^p(v)(NF^v)y\l + o(l)), for u G (-E^, 0). 

For u > 0, ECrtAr(u) = 2ECrtAr(0)(l + o(l)), obviously. Inserting back the definitions 
(7.1) completes the proof of Theorem 2.17[ c,d). 



7.2. Sharp asymptotics at and beyond the edge. We first rewrite the one-point 
correlation function p^ using the Hermite functions j G N, given by 

cj> j (x) = (yj\Vn)- 1 / 2 H j (x)e-£, 



(7.4) 

where Hj, j € N are Hermite polynomials, Hq(x) = 1 and Hj(x) = e x (~^y e ~ x ■ The 
Hermite functions are orthonormal functions in K with respect to Lebesgue measure. From 
|Meh91| . pp. 128 and 135, it follows that 



p N (x) = N~ 1 / 2 (S N (x) + a N (x)), 



(7.5) 
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where 



S n (x)=Y / $(VNx) + (-) c^N-ii^Nx) e(VNx-t)<f> N (t)dt, 
i=o J -°° 



e(x 



a N (x 



if N = 2m + 1, 
if N is even. 



(7.6) 
(7.7) 
(7.8) 



sign(x)/2, 

Ifomi^Nx^J^fomWdt}- 1 , 

The factor A^ 1 / 2 in (frTsl) comes from a change of variables and the fact that the one-point 
correlation function R± in |Mch91| is not normalized to be a probability density. 

Using the Chistoffel-Darboux formula QMch91], p. 420), the first term of (7.6) satisfies 
JV-l 

tfi^Nx) = N4> 2 N (VNx) - y/N(N + IWn.^VNxWn+^VNx), (7.9) 

i=0 

so that pn(x) depends on 4>n-i, 4>n and 4>n+i only. We now state the Plancherel-Rotach 
asymptotics of the Hermite functions in domains of our interest as a lemma since it will 
be our main tool from now on. We use tp, and h to denote the functions 

ij,{x) = \ x 2 - 2| 1/2 , Ti(v,-=)(x) 



V2 J 



sft 



h(x) 



X-V2 


1/4 


x + s/2 


1/4 


x + y/2 


+ 


x-V2 





(7.10) 



Lemma 7.1. There exists a 5q > 0, such that for all < 5 < 5q the following holds 
uniformly for x in the given domains. 

(a) Fory/2-5 <x<yj2 + 8, 



4> N (VNx) 



1 

: (2A0V 4 
x 



x + V2 
x-V2 



1/4 



V2 



+ V2 



1/4 



1 



|/ w (x)| 1 / 4 Ai(/^(x))(l + 0(iV- 1 )) 



■Ai'(/ iV ( a ;))(l + 0(iV- 1 )) 



(7.11) 



\fN(x)\V* 

where /at(x) = A r2//3 {|/i(t>, ^)(x)} 2 ^ 3 , and Ai(x) is the Airy function of first 
kind, Ai(x) = | cos + fcc)dt. 



(b) For x > y/2 + 5, 



4>nWnx) = , * w (l + o^- 1 )). 



(7.12) 



\/47r^/2iV 

Proof. The lemma follows from [D KM + 9"9"1 IDG07j . Formulas closest to our formulation 
can be found in |DG07| . pp. 20-22. Under their notation, our c ase has Cjy = y/2N, 
h N (x) = 4 and m = 1. For definitions of / see (2.15), (2.16) of [DKM+99) . For the 
original Plancherel-Rotach asymptotics for the Hermite polynomials, the reader can also 
check |PR29| or |Sze81j . □ 



For the rest of this section, we suppose that 5 is small enough so that Lemma |7.1| 
holds. We write ~ bjy, if cln = &Af(l + o(l)) as N — > 00. For sequences ajy(S), &tv(<5) 

which depend on 8, we write ~ if lim^o lini/v-»oo ^{S) /bjsf(S) = 1. Finally, if 
a^v < C(5)on for some C(<5) < 00, we write = Og(b^f). 



We first analyse the integrals appearing in (7.6) and (7.8). We define 

/*oo /*oo 

J N {x) := / e{^Nx-t)<t) N (t)dt = N 112 I e(x - s)<p N (sVN) ds. 



(7.13) 
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(f> N (x) dx = 2 <j> N {x) dx ~ 2(2N)~ 1/4 . 
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(7.14) 



(b) If x > \/2 and N — > oo over odd integers, then Jn(x) = 0(N c e NIl ^"'V2^ x ^ M 

(c) If x > y/2 and N — > oo over even integers, then Jn(x) ~ (2AQ -1 / 4 . 

(d) Ze£ a^v > be such that limAr_-. 00 a^N 2 / 3 = 0. T/ien 



Jjvfv 2 + a at) ~ 



i(2JV)-V* 
-i(2iV)- 1 /4, 



as iV —7- oo over even integers, 
as N — )■ oo over odd integers. 



(7.15) 



Proof. Claim (a) can be derived directly from Proposition 4.3 of [D G07| . 
The Hermite functions are even (odd) for iV even (odd). Therefore, 



J N (x) 



even, 



From Lemma 



AT 1 / 2 sign x Jj^^(sv^V)ds, iV odd. 
7.1 b) it follows that for x > \[2 and some c, C < oo 

^>75)M ds = 0(7V c e- iV/l( "'75 )W ). 



</> N (sVN) ds < CiV c 



(7.16) 



(7.17) 



To prove the last equality we used the fact that h(v, 4^) is strictly increasing on (y2, oo) 
and Laplace's method. Claims (b), (c) are then direct consequences of (7.17) and claim (a). 



In view of claim (a), (7.16) and (7.17), to prove (d) it suffices to show 

8 1 



V2+S 

^(sVN^s ~ ^2- 1 / 4 iV- 3 / 4 . 



(7.18) 



To this end we use Lemma |7.1[ a). We first linearise the function /jy appearing there. It can 
be proved easily from its definition that for all N > 1, uniformly over x G {y2 — S, y2 + S) 

f N (x)l2^N^(x-V2-), 
f' N {x) I 2 l ' 2 N 2 /\ 
Hence, after substitution /tv(x) = z we obtain using Lemma |7.1[ a), 



(7.19) 



. ,/ ]v (v / 2+<5) , 

/ ||2 3 / 2 |V4( JV 2/3 2 l/2 ) l/4 Ai( , 



V2+S l 

(P N (sVN)ds ~ T-TTTJm , . 



(7.20) 



- |2 3 / 2 |- 1 /4 (iV 2/3 2 l/2 ) -l/4 Ai / ( ^ } j 2 -l/ 2iV - 2 /3 

The limits of the integral converge to and oo respectively. Using the well-known identity 
J °° Ai(z) dz = |, the claim (7.18) follows. This completes the proof of the lemma. □ 

We can now compute ECrtjv(iO for u < —E^. By (7.2), (7.5), (7.9), 

4 „™ 4 



e iV ^)r,(x)dx =: CW.pJV" 1 /* (7.21) 



where 



Ti(x) = «at(x), 
T 2 (x) = N(f> 2 N (VNx), 

T 3 (x) = —\J N(N + IWn^VNxWn+^VNx), 

/ATnI/2 / _ 
T 40*0 = IT <j> N -i(VNx)J N (x). 



(7.22) 
(7.23) 
(7.24) 

(7.25) 
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Proof of Theorem 2.11(a). We first estimate the four integrals Ii(v), . . . , I&(v) for v > y/2. 



Lemma 7.3. For N even I\(v) = 0. When N — > oo over odd integers, then 

h(v,^)(v)-^(v) 



h(v) 



-AT(J!>(«) + / 1 («,-^)(«)) 



(7.26) 



claim we fix o" < v — \/2 and use Lemmas 7.1 ^b), 7.2 a). Then, 



Proof. The first claim follows directly fro m the de finit ion ( 7.8 ) of a at. To prove the second 

(7.27) 



h{v) 



(2JV) 1 /* p /i( x ^^ (a) ~ (N_1)/l( "'75 )(a!+) 



dx, 



or 



(4vr) 1 /2( 2 (iV- 1))V4 

where x + = x^N/(N - 1) ~ x(l + l/(2iV)). Expanding h{v, ^)(x+) ~ ^)(x) + 
xl\(v, ■^)'(x)/(2N), we obtain the lemma by Laplace's method. □ 

The integrals l2(v) , ^3(^) are negligible for v > \/2: 

Lemma 7.4. There exists a c < 00 suc/i f/iai Ij(t>) = 0(N c e~ N( '~' t ' ( ' V '' +2Il( - v, ^ ( - v ^) /< 
i = 2,3. 

Proof. Observe that T2, T3 contain a product of two Hermite functions. By Lemma [7.1[ b), 
these behave like 0(N c e' Nh{v '^ ){x) ) for x > a/2, which implies the claim. □ 

Lemma 7.5. (a) // N — > 00 over odd integers, then 
for some c < 00. 



(b) If N — > 00 over even integers, then Ia(v) behaves like the right-hand side of (7.26). 



Proof. Claim (a) follows from Lemma |7.2[ b) and the same reasoning as in the previous 
proof. For claim (b), we find using Lemma |7.2|(c ) and Lemma |7.1[ b) that I4 is asymptot- 
ically equivalent to the right-hand side of (7.27) and the claim follows. □ 

Theorem 2.17[ a) follows directly from the previous three lemmas and definitions (7.1), 
(7.10). Observe that the dominant contribution comes from I\ for N odd and from I4 for 
N even. □ 

Proof of Theorem 2.11(b). We need now compute h(v), . . . ,Ia(v) for v = y/2. We split 
all these integrals into two parts: over (y/2, y/2 + 8) and over (v2 + 8, 00). For the second 
interval we can use Lemmas 7.3-7.5 We need thus analyse the integrals over the first 

interval, if = J^ +S Ti(x) dx, i = 1, . . . , 4. 



Lemma 7.6. For N even h{y/2) = 0. When N — > 00 over odd integers, then 



6 T s S 



Ai(0) 



,N<f>(V2) 



21/2^5/6^/(^/2) 



(7.28) 



Proof. Due to Le mma 7.3 it suffices to prove the second '~' relation only. By Lemma 7.1 a) 
a), setting 



and Lemma 



7.2 



T 5 5 1 



\/2+5 



v/2 



N4>{x) 



X+ + y/2 



1/4 



x+- y/2 
X+-V2 
x+ + yf2 



\f N (x+)\ 1 / 4 Ai(f N (x + )) 
V4 |/iv(x + )r 1 / 4 Ai'(/ w (x + ))Ux 



(7.29) 
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If it were not for the fpj(x+) in the argument of the Airy function, this integral could 
be analysed trivially by Laplace's method, the main contribution coming from neighbour- 
hoods of size O(N^) of y/2. However, due to fTl9| ), f N (x+) converges to uniformly 



over such neighbourhoods. Hence, by a simple extension of Laplace's method, using (7.19) 
several times, 

S £ Al (0) MV2) (J on\ 

1 21/2^5/6^(^2) 6 • V' M > 

This completes the proof. □ 
Lemma 7.7. As N — > oo, 

I 2 (V2) = 0{N- 1 / & e N ^) and I 3 {V2) = -I 2 (V2)(1 + C^iV" 1 )). (7.31) 



Proof. As before, /2(\/2) ~ if- To estimate /| we use Lemma 7.1 a) again. The same 
computation as in the proof of the previous lemma gives 

sl Ai 2 (0)V2_ ^ (V5) 

2 N^F^V2) ^ ' 



To prove the second claim, it suffices to observe that on (y/2, y/2 + 5), by Lemma 7.1 'a), 

^ Xy/N) = = 1 + 0(N^). (7.33) 

□ 

Lemma 7.8. (a) When N — > oo over even integers, then 

h(V2) I 4 I \g n ( P ) (7.34) 



where Gn{p) denotes the right-hand side of (7.28). 
(b) When N — > oo over odd integers, then 

5 T s 6 1 



I 4 (V2)l4^--G N {p). (7.35) 



Proof. The lemma follows by the same computations in Lemma 7.6 Again, the dominant 
contribution comes from neighbourhoods of size OiN- 1 ) of v^T On such neighbour- 
hoods, using Lemma |7.2|^d), the integral Jjy(x) appearing in T±{x) can be approximated 



by cN x / 4 , where c depends on iV being even or odd. □ 



Theorem 2.17[ b) follows now directly from the previous three lemmas and definitions 



(7.1), (7.10). Observe that the dominant contribution comes from I\ and I4, since I 2 and 
I3 mutually cancel. Moreover, the combined contributions of I\ and I4 do not depend on 
N being odd or even. □ 

Proof of Corollary \2.18[ By Theorem |2.5[ for all e > 0, 

ECrtjv(it) > ECrtjv.ofa) > ECrtjv(u) - AHE Crt at,i (u) 
> E Crtjv(tO - Ne N{ei ^ u)+£) . 

The corollary then follows from the fact that for u < — E^, @q iP (u) > &i p (u). □ 
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Appendix A. Large deviations for the largest eigenvalues of the GOE 

In this appendix we extend Theorem 6.2 of |BDG0l| . proving a LDP for the k-th. largest 
eigenvalue of the Gaussian Orthogonal Ensemble. This results might be of independent 
interest. Its proof follows the lines of [BDG01]. 

Let X = Xn be a N x N real symmetric random matrix whose entries X^ are inde- 
pendent (up to the symmetry) centered Gaussian random variables with the variance 

EXj = a 2 N~ 1 {l + 8 ij ), (A.l) 

and let Ai < • • • < An be the ordered eigenvalues of X. Note, that in this appendix we 
use the notation that is more usual in the random matrix theory, that is the eigenvalues 
are numbered from 1 to N, not from to A^ — 1 as in the rest of the paper. We show the 
following LDP. 

Theorem A.l. For each fixed k > 1, the k-th largest eigenvalue \N-k+i °f ' X satisfies a 
LDP with speed N and a good rate function 

I oo, otherwise. 

Proof. We first recall some know fact about the distribution of Aj's and introduce some 
notation. The joint law of Aj, 1 < i < N, is given by 

N N 

Q N (d\i, . . . ,d\ N ) = Z N (aY 1 \\ - Xj\ JJexp ( - Af^dAjlAi^-^Ajv ( A -3) 

l<i<j<N i=l 

The distribution of unordered eigenvalues Qn is given by the same formula without the 
final indicator function and with Z^(a) replaced by Z^(a) = N\Z^(a). By Wigner's 
theorem, the spectral measure Ljy = A^ -1 YliLi ^\ °f -^n converges weakly in probability 
to the semi-circle distribution 

p(dx) = (27ra 2 )- V(2^) 2 - x 2 l\ x \<2* dx. (A.4) 

For AcKwe denote by V(A) the space of all Borel probability measures on A endowed 
with the weak topology and a compatible metric d. By Theorem 1.1 of [BG97], the spectral 
measure Ln satisfies a LDP on V(M) with the speed A^ 2 and a good rate function whose 
unique minimiser is the semi-circle distribution (A.4). 

We now start proving Theorem |A.l If.(x;a) is obviously good rate function. Since 
Ik(x,a) is scale invariant, that is 

I k (x;a) = h{x/a; 1), (A. 5) 

we can assume that a = 1 as in [BDG01], and omit a from the notation. Note that this 
value differs from the value a = 2 -1 / 2 used in the rest of the present paper (cf.( |2.6| )). In 
particular, for a = 1 (cf. ( 3.17| )) 

^ = ^(l) = (2V2)~(f)-~nr(l + i). (A.6) 

1=1 



To show Theorem A.l it is sufficient to prove 

limsup — log Qn (XN-k+i < x) = — oo for all x < 2, (A. 7) 
and, since Ik(x; 1) is continuous and strictly increasing on [2, oo), 

limsup — log Q n (XN-k+i > x) = —Ik(x;l) for all x > 2. (A. 8) 



We first prove (A. 7). Suppose that Xjy^k+i < x for some x < 2. Then Ln((x, 2]) < 
(k — \)N~ l . Since p((x,2]) > 0, for N large enough there exists a closed set A C V(M.) 



RANDOM MATRICES AND COMPLEXITY OF SPIN GLASSES 24 

such that p £ A and {^N-k+i < s} C i. By the LDP for the spectral measure Ljv, 
Qn(A) < e~ cN2 for some c > 0, which concludes the proof of (A.7). 
We now prove the upper bound for (A.8). By Lemma 6.3 of [BD GOlj . 



Q N f max | Ai| > AfJ < e~ NM /y for all M large enough and all N. (A. 9) 

Writing 

/AT n / N 

QN{^N-k+i > x) < QN{max\\i\ > M) + Qwy^N-k+i > £,max|Ai| < M), (A. 10) 
the upper bound follows easily provided we show that for all M > x > 2 

limsup ^- logQAr(niax |Aj| < M, Xx-k+i > x) < -Ik{x;l). ( A -H) 

AT->oc ^ 1=1 



To show (A.ll) we introduce some additional notation. Let Qjy-fc ^ e a measure ° n 
R N ~ k given by 

0jv_fc(A G •) = Qiyr_*((l - kN'^X E •)• (A.12) 

We set 

Cn -{ 1 -n) (A - 13) 

and for i£l and \i G V(M) we define 

f z 2 
$(z,n) = / log\z - y\n(dy) - — . (A.14) 

It was shown in [BDG01, p. 50] that §(z,p) is upper semi-continuous on [— M, M] x 
V([-M, M\) and continuous on [x, y] x Pfl-M, M]) for all M, x, y G M. such that y > x > 
M > 2. 

Using ( |A.3 ) and this notation, we can write 
QAr(max|Ai| < M, Xjsf-k+i > x) 

N N—k 
J[x,M] k i= N k+ l J[-M,M] N - k t=l 



X II l^i — ^j|lAi<--<A ]v 

l<i<j<N 

- ckN W^ky.L M]k II \x i -x 3 \d\ N . k+1 ...dx N 

x / et N - k ^?= N - k+1 H\,LN- k )QN_ k{dXl ,...,dX 

J[-M,M] N ~ k 



(A.15) 



where the factor N\/(N—k)\ comes from replacing Qn by Qn. Let B(p, 5) denote the open 
ball mV(R) of radius 5 > and center p. We write B M {p,5) = B(p,5)C\V{[-M,M}). On 
the domain of the integration \Xj -Xj\ < 2M and e (N-k)*(\i,L N - k ) < (2M) N ~ k . Therefore 
by splitting the second integral, (A.15) is bounded from above by 



G 



k N\ - - - ;/ 

N~. 



f2M ^(fc-l)/2f/ / (W-fc)BU P((6BMM *M d \ 

( Ar - fc ) ! U A / (A.16) 



+ 



(2M) Ar - fc A >- fc (^-^ 5(^,5))}. 
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To control the measure Q^_ k , observe that for all functions h : K - 
norm at most 1 and N >2k, 
N-k 

(N - k)- 1 {h{{l ~ kN-^Xi) - h(Xi)} 



^ »t-1 N-k , , 
< cJS max A,- 



i=k 
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of Lipschitz 
(A.17) 



for some c E (0,oo) independent of N and k. It follows from (A. 9) that the laws of 
LN-k under Q^-k and Q%_ k are exponentially equivalent as N — > oo. Therefore, by 
Theorem 4.2.13 of |DZ98| . L]^_ k under Q%_ k satisfies the same LDP as L^_ k under 
Qx-k- Hence, the second term in (A. 16) is exponentially negligible for any 5 > and 
M < oo. This implies that 



lim sup — log Qn( max | \ \ < M, A N-k+i > x) 



< lim sup — log C% + k lim sup $(z,/i). 

Y->oo N 810 ze [ x> M],neB M ( P ,S) 



(A.18) 



The same reasoning as on p. 50 of [BDG01] implies that the second term in (A.18) equals 
-k(l/2 + h(x; 1)). From definition flA.13[ ) of C% and from ( [AT6| ) , it is easy to obtain 
limjv-s.oo N~ l log C^r = k/2. Combining these two claims, we can bound the left-hand side 
of ( [A. lip by -k(l/2 + I 1 (x; l)) + k/2 = kh(x; 1). This completes the proof of ( |A.11[ ) and 
thus of the upper bound for (A.8). 

To prove the complementary lower bound wefixy>x>r>2 and 5 > 0. By a similar 
computation as in ( |A.15 ) we obtain 

QN^N-k+i >x)> Qn(^n G [x,y],.. -,XN-k+i e [x,y],max|A;| < r) 

i=i 



N 



N 



n e -^?/ 4 dA 4 n 

i=N-k+l N-k<i<j<N 



Xi — Xj 



(A.19) 



,(N-k) *Q*,L»-k)QN_ k ( dXl , dX N _ k ) 



>KC k N e W (k(N-k) inf ^z,vi))Q%_ k {L N - k € B r (j>,8)), 

for some = x, y) > 0. 

Using the LDP for the measure L^_ k under Q N -_ k we see that 



Q%_ k (L N _ k £ B r (a,S)) *t22> q. 



By symmetry of Qn(') and by the upper bound in (A.8), 



Q%_ k (L N ^ k i P((-r,r))) < 2Q N . k (X N ^ k > r) 



Therefore, using the behavior of again, 

1 k 
limM — logQ N (X N _ k+ i> x) > - + k inf ®{z,(j). 

A->oo Jy 2 zG[x,y},ueB r (p,6) 

Letting now 5 — > and then y \, x, using the continuity of $>(z,fj,) in the used range of 
the parameters, we obtain the desired lower bound. □ 



(A.20) 
(A.21) 

(A.22) 
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